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Abstract. We use a Lucky Imaging system to obtain I-band images witii improved angular resolution on a 2.5m telescope. We 
present results from a 10-night assessment campaign on the 2.56m Nordic Optical Telescope and quantify the performance of 
our system in seeings better than 1.0". In good seeing we have acquired near diffraction-limited images; in poorer seeing the 
angular resolution has been routinely improved by factors of 2.5-4. The system can use guide stars as faint as 1=16 with full 
performance and its useful field of view is consistently larger than 40" diameter. The technique shows promise for a number of 
science programmes, both galactic (eg. binary candidates, brown dwarfs, globular cluster cores) and extragalactic (eg. quasar 
host galaxies, damped Lyman-a absorbers). 
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1. Introduction 

We have demonstrated in recent work the recovery of essen- 
tially diffraction limited I-band images with a novel imag- 
ing system on the 2 . 56m Nordic Optical Telescope (NOT) 
fealdwin et al."2001; 'Tubbs et alj |2002l l2003t iMackav et all 
t2004b). Our Lucky Imaging system (LuckyCam) takes a se- 
quence of images at >10 frames per second using a very low 
noise L3CCD based conventional camera. Each short expo- 
sure suffers different atmospheric turbulence effects, resulting 
in rapid variations in image quality. To construct a final image 
we select and co-add only those frames which meet a quality 
criterion. By varying the criterion we can trade off sensitivity 
against ultimate resolution. 

Lucky Imaging is a passive technique, so useful data is 
taken as soon as the telescope is pointed correctly. LuckyCam 
may be used in the same manner as any normal CCD camera 
system, with no special setup or operational requirements. 

Lucky Imaging shares the goal of Adaptive Optics (AO) 
systems - to enhance ground-based telescope resolution. AO 
systems are currently in use at several telescopes, and give 
very valuable resolution in the infrared. However, we are not 
aware of any systems routinely giving excellent resolution over 
a large (>5") field of view at wavelengths shorter than 1/im. 

In this paper we show that our prototype Lucky Imaging 
camera can give a very valuable increase in I-band resolution 



* Based on observations made with the Nordic Optical Telescope, 
operated on the island of La Palma jointly by Denmark, Finland, 
Iceland, Norway, and Sweden, in the Spanish Observatorio del Roque 
de los Muchachos of the Institute de Astrofisica de Canarias. 



over a large area with only modest signal loss due to frame se- 
lection. Using data taken over 10 nights in 2000-2004 we ex- 
plore the system's performance in I-band on a 2.5m telescope. 
This investigates a small part of the possible Lucky Imaging pa- 
rameter space, which can be defined in terms of camera frame 
rate, passband and telescope size. 

In most observations presented here the camera is run ap- 
proximately 3x too slowly to fully sample changes in the at- 
mospheric turbulence and so we cannot demonstrate properly 
the full potential of Lucky Imaging. However, the very signifi- 
cant increase in resolution available with the current system is 
valuable for many science programmes. 

In section |2l we describe our camera and data reduction 
techniques. In section |3l we describe the point spread func- 
tion improvements produced and the specific data reduction 
techniques we use. In section Owe detail the results of Lucky 
Imaging trials to date. Section |5] addresses the effect on limit- 
ing magnitudes of the Lucky Imaging frame selection; section 
|6ldiscusses the guide star and hardware requirements. 



2. The Lucky Imaging Technique 

The Cambridge Lucky Imaging system (LuckyCam) is based 
on an E2V Technologies L3CCD read out with a 4 MHz pixel 
rate and mounted at the focus of a simple reimaging camera. 
The on-chip gain stage of the L3CCD raises the signal from 
incoming light sufficiently to allow individual photons to be 
detected with good signal to noise, even at hi g h frame rates 
- for more details see e.g. M ackav et a l] |l2003l) : lBasden etalJ 
I2OOI) : Mackav et aL(.2004ai) . 
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Fig. 1. Lucky Imaging of the core of M15. The left panel is the output from an autoguider system (FWHM 0.63"); the right shows 
the 10% frame selection Lucky output from the same exposures (FWHM 0.26"). Both images are on the same grayscale. The 
increased resolution makes more faint objects visible in the Lucky image. 



Typically the camera is run in full-frame (552x512 pixel) 
12 frames per second (FPS) mode, although the frame size 
may be reduced for higher frame rates. The image scale is 
0.04"/pixel, giving a field of view of 22x20.5 arcseconds. 
This slightly undersamples the 0.08" FWHM I-band diffrac- 
tion limited PSF that would be produced by the Nordic Optical 
Telescope. 

The frame selection algorithm, implemented (currently) as 
a post-processing step, is summarised below: 

1 . A Point Spread Function (PSF) guide star is selected as a 
reference to the turbulence induced blurring of each frame. 

2. The guide star image in each frame is sinc-resampled by a 
factor of 4 to give a sub-pixel estimate of the position of the 
brightest speckle. 

3 . A quality factor (currently the fraction of light concentrated 
in the brightest pixel of the PSF) is calculated for each 
frame. 

4. A fraction of the frames are then selected according to their 
quality factors. The fraction is chosen to optimise the trade- 
off between the resolution and the target signal-to-noise ra- 
tio required. 

5. The selected frames are shifted-and-added to align their 
brightest speckle positions. 

Faint guide star PSFs are affected by photon shot noise, 
which can lead to an image being positioned on the basis of 
a noise spike rather than an actual bright speckle. This noise 
is reduced by convolving the faint reference star image with a 
theoretical diffraction-limited PSF. 

An imp lementation of the Drizzle algorithm 
jpruchter & Hook ■2OO2I) is used for the image alignment 



step. The algorithm resamples the images by a factor of two 
and minimises the information lost by our slightly under- 
sampled pixel scale. Drizzle is especially suitable as we are 
summing many frames displaced by random non-integer pixel 
offsets. 

Figure Ogives a typical example of the improvements ob- 
tained with Lucky Imaging. 

3. The PSFs produced by Lucky Imaging 

Examples of the general form of the Lucky Imaging PSF are 
given in fig.|3 As with adaptive optics images, the radial shape 
takes the expected form for an ima ge with partial ly compen- 
sated Kolmogorov turbulence (e.g.. iHardvl lll998l) and refer- 
ences therein) - a wide halo and a central compact core. As 
the selection of frames is made more stringent the fraction of 
light in the compact central core progressively increases. 

The Strehl ratio (commonly used as a high-resolution imag- 
ing performance measurement) is the peak value of a PSF 
divided by the theoretical diffraction-limited value. Figure |2 
clearly shows that Lucky imaging offers a substantial Strehl 
ratio improvement over both autoguider and shift-and-add (ef- 
fectively 100% selection) systems. We have routinely measured 
I-band Strehl ratios in the range of 0. 15-0.2 at high frame rates 
in good seeing on the 2.5m NOT. 

We note that the NOT is not designed or caUbrated to pro- 
duce diffraction limited images. Although the random phase 
variations in the atmosphere can compensate for telescope mir- 
ror (or focusing) errors they do so at a lower probability than 
the production of a good image with a perfect telescope. In the 
case of the NOT small scale mirror irregularities limit the peak 
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Fig. 2. Average radial Lucky Imaging profiles. Data for the left figure was taken at 30Hz in 0.6" seeing; the right figure shows the 
results of 12Hz imaging in 0.55" seeing. 100% selection is simple shift-and-add. When selecting 1% of frames the Strehl ratio is 
more than doubled relative to the 100% selection, while the light from the star is concentrated into an area approximately four 
times smaller. 
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Fig. 3. FWHM of target stars 1-3" from the reference star. Lines correspond to resolution increase factors of Ix, 2x, 3x & 4x. 
Most of the scatter in the FWHMs is caused by the range of target star to guide star distances. Vertical crosses are runs taken at 
12Hz in 2003 & 2004; diagonal crosses are runs taken at 18Hz in 2001 & 2002 and circles are runs taken at 36Hz in 2003. 



Strehl in the absence of an atmosphere to around 0.2. Optics 
designed for high Strehl ratios would increase the probabil ity 
of obtaining a high quality frame substantially llTubbJ2004h . 

We cannot directly compare the Lucky Imaging perfor- 
mance to a standard imaging system as we are not able to take 
images simultaneously with the two systems. Rapidly changing 
seeing can then lead to a bias. To avoid the bias in this paper 
we apply a simulation of the NOT autoguider system to our 
LuckyCam data, effectively giving a simultaneous test of the 
two systems. 

The quoted seeings are measured from the FWHM of 2D 
Moffat profile fits to simulated autoguider images and agree 
with those measured at intervals by conventional means during 
the observations. They are thus measured at an effective wave- 
length of ~ 850nm and should be increased by approximately 
10% to convert to a standard seeing quoted at 550nm. All other 



FWHMs detailed in this paper are similarly measured from 2D 
Moffat profile fits. 

4. Lucky Imaging Performance 

The present results are based on observations taken during 
10 nights on five observing runs in May, June & July, 2000- 
2004. The 2.56m Nordic Optical Telescope on La Palma in 
the Canary Islands was used for all observations, which were 
principally made in several different fields in the cores of the 
Globular Clusters M3, M13 and M15. The 2001 -2004 dataset is 
split into 42 short (2-3 minute) runs, totalling ~ 1 32,000 frames. 
All fields were observed in I-band. 

We here investigate science-target performance - i.e. that 
for a star near the guide star In each of the fields we chose 
a target star for PSF measurement at 1-3" separation from the 
guide star. 
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Fig. 4. The factors by which FWHM and FWHEF (full width at half enclosed flux) are improved at slow frame rates (12FPS), 
at three seeings. Obtained from linear fits to figure|5]and similarly derived results for the FWHEF. Starting at the leftmost point 
of each fine, datapoints correspond to selecting frames at the 1%, 2%, 5%, 10%, 20%, 30%, 40%, 50%, 75% and 100% levels 
respectively. 




The improvement in FWHM obtainable in different atmo- 
spheric conditions and with diff'erent percentage selections is 
shown in fig.|3] Much of the scatter in values is due to the range 
of distances of target stars from the reference star (section lTTl . 

Figure |3] shows that under a wide range of conditions the 
resolution is improved by factors as large as x4 in the most 
stringent selections in all of the 42 observations reduced. Less 
stringent selections give smaller improvement. There is an 
approximately linear correspondence between the autoguided 
long exposure seeing and the resolution attainable. This sug- 
gests that, at least over the ~200 second timescale of these 
observations, the standard measures of seeing are a reason- 
able guide to the atmospheric turbulence statistics. We can thus 
adopt the seeing as our standard measure of the atmospheric 
turbulence strength in a particular run. 

Figurel^shows the effects of selecting dififering fractions of 
images using empirical fits to the 12FPS data. Although limited 
by the slow frame rate, in 0.5" seeing the most stringently se- 
lected images are within ~0.06" of the 0.08" diffraction limit 
and Strehl ratios are >0. 1 . 

The advantages of the full-frame field of view led us to ob- 
tain most of our data in LuckyCam's relatively slow 12 or 18 
FPS modes. To obtain images with light reliably concentrated 
into bright single speckles, unblurred by image motion, we 
must oversample the atmospheric coherence time. This would 
typically require > 40 FPS at the NOT. However, the image 
quality varies on every timescale slower than the coherence 
time and so image selection can always be expected to improve 
the resolution - even if we use very slow frame rates. 

To acquire a high-quality frame we require both an excel- 
lent point spread function and a stable atmosphere for longer 
than our frame integration time; this occurs with increasing 
probability as we increase the frame rate. As expected, the 
faster frame rates presented in fig.|3lon average give a greater 
resolution improvement than the slower rates. In good seeing 
the output resolution at fast frame rate is within 0.03" of the 



Fig. 5. A 1 % selected image of ^ Bootis with a Strehl ratio of 
0.26 and a FWHM of ~0. 1", taken in 0.42" seeing. 

diffraction limit of the telescope, limited at least in part by 
LuckyCam's slight PSF undersampling. 

Although the resolution can be adjusted by altering the 
number of frames selected, in poorer seeing the probability 
of a superb frame is greatly reduced. Empirically this appears 
to limit the LuckyCam resolution increase to a factor of ~3-4 
when the seeing is poor 

The FWHM gives a good estimate of the system's perfor- 
mance for resolution enhancement - i.e. the separation of two 
closely separated objects. However, for many observations (no- 
tably precision photometry) the degree of light concentration 
within the larger halo is also important. Figure |4] also details 
the Full Width at Half Enclosed Flux (FWHEF) performance. 
The FWHEF is the aperture size which contains half the light of 
the star and thus is a better measure for crowded observations 
where photometry is to be performed. The current LuckyCam 
improves the FWHEF by up to a factor of two, corresponding 
to concentrating half the light from a star into an area four times 
smaller than in seeing-limited images. 

With a faster camera higher resolutions have been reached. 
Figure |5] shows an I-band image of the 0.8" binary ^ Bootis 
with a Strehl ratio of 0.26, taken in 0.42" seeing. The first 
Airy ring is clearly visible at a radius of ~0.1", showing the 
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Fig. 6. An example set of isoplanatic patch extrapolations. 
Measured Strehl ratios of a number of stars (from a 0.5" seeing 
run at 1 % selection) are plotted as a function of distance from 
the reference star; three empirical models are fit. 



point spread function is indeed diffraction-limited in width. 
This dataset was taken at 200Hz, approximately 5x faster than 
the coherence time of the atmosphere requires. 

4.1. Isoplanatic Patch 

Resolution enhancement systems such as LuckyCam provide 
useful turbulence correction only within some angular distance 
to the guide star This distance is commonly quantified by the 
isoplanatic patch radius - the radius at which the Strehl ratio 
of stars is reduced by a factor of e ' from those very close to 
the guide star We here use this definition of our useful field 
size, although noting our low frame rate wide-field data has 
generally smaller Strehl ratios (of order 0.1) than are usually 
used for this measurement. 

The difference in Strehl ratio between stars close to the 
guide star and those up to 25" off-axis (our field limit) can 
be measured for several of the globular cluster fields in our 
dataset. The isoplanatic patch radius is then calculated from the 
point at which an empirical function fitted to the Strehl ratios 
drops to e"' of its on-axis value. The models (fig.|6j suggest 
that in five fields of Ml 5 taken over two nights in 2003 the iso- 
planatic patch size ranged from 17-30" in radius. Selecting a 
larger fraction of the frames gives an improved patch size (up 
to 10" wider) as the degree of turbulence correction achieved 
on-axis is decreased along with the output resolution. 

It is clear that the atmospheric turbulence corrections pro- 
duced even with a low frame rate camera have a remarkably 
large effective radius, in most cases larger than LuckyCam's 
field of view. These results (from data taken in 2003) agree with 
earlier Lucky Imaging results presented in Tubbs et al. (2002), 
suggesting that the large patch size may be inherent to the tech- 
nique over a wide range of conditions. 

The isoplanatic patch for speckle imaging is expected to 
be pro portional to rn, the standard atmospheric coherence 
length llVernin & Munoz- Tunonllll994 l) and references therein), 
ro is known to vary on relatively short timescales (fig. Q. If 
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Fig. 7. An example of fast seeing variations. These images were 
taken with an 83ms exposure time in the core of M3 in average 
0.85" seeing - 0.33 seconds apart. This is not a results of the 
short timescale statistical fluctuations in image resolution ex- 
pected for a constant ro - the 5-second averaged seeing changed 
by approximately a factor of two between these two images. 
Although extreme, these fast variations are present in most of 
the runs analysed here. Each bright pixel in these images is a 
detection of at least one photon. 
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Fig. 8. The SNR (both measured and modelled) of an I=+19m 
point source in a 16 minute exposure as a function of the Lucky 
Imaging improvement in FWHM. The aperture size is set to 
the PSF diameter at half enclosed flux. Starting at 100% at the 
upper end of each line the data points represent frame selection 
fractions of 100%, 75%, 50%, 40%, 30%, 20%, 10%, 5%, 2% 
and 1%. Note that these observations were made in the central 
core of M15 and have ~ 4x higher background than an empty 
field - i.e. the SNRs are reduced by a factor of two. 



LuckyCam selects periods when ro is larger (and not just pe- 
riods when the turbulence induced phase errors randomly sum 
to a small value) it would be expected that those periods would 
have larger isoplanatic patches, as we observe. 

5. Signal-to-noise Ratio (SNR) Considerations 

An improvement in resolution allows smaller photometric aper- 
tures to be used. A smaller aperture contains fewer pixels and 
so less sky noise and (here negligible) detector noise, giving an 
increased SNR. 

In fig. |8l we compare modelled LuckyCam SNR perfor- 
mance with optimised apertures (based the performance de- 
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Fig. 9. The factor by which the target star FWHM is improved 
as a function of the guide star flux at a 10% selection. The 
dashed lines show results based on simulating a faint star by 
adding photon & L3CCD multiplication noise to a bright star's 
PSF; also shown (with crosses) are measurements using actual 
stars. 

tailed in section 0} to measured results for an I=+19m star 
in the core of M15. As for this star the noise in these images 
is background dominated, the SNR is measured by comparing 
the signal contained inside an aperture to the RMS noise in a 
nearby empty (background) aperture. 

The 0.55" seeing model agrees well with the measured 
SNR in high-resolution observations but diverges by ~25% at 
lower resolutions. This is due to increased crowding (our fields 
are in the cores of globular clusters) leading to an increased 
background. This aside, the standard CCD photometric noise 
estimates model LuckyCam's performance well, once the ap- 
plicable aperture size changes are taken into account. 

Although a 1 % selection of frames would be expected to 
decrease the SNR by a factor of VlOO = 10 for a star in an 
empty field, the smaller aperture sizes that can be used reduce 
the decrease to ~ 6.5x in both seeings in fig.|8]. At the 50-75% 
selection level the SNR matches that available in a standard 
autoguided image, with resolution increased by a factor of ~ 2. 

We emphasise that Lucky Imaging 's ability to increase res- 
olution by large factors even in poor seeing eff'ectively in- 
creases the usable observation time at the telescope. For ex- 
ample, if a particular observation requires 0.5" seeing selecting 
a few percent of frames in 1.5" seeing will produce usable data 
in telescope time that would otherwise have been lost. 

6. Lucky Imaging Requirements 

In this section we detail the two main requirements for a suc- 
cessful Lucky Imaging observation - a sufficiently bright guide 
star and a sufficiently fast camera. 

6.1. Reference Star Flux 

We simulate the effect of using a faint guide star by rescaling 
a bright (> 2000 photons / frame) star PSF to a specified ffux. 



taking Poisson and L3CCD multiplication noise into account. 
As a check of the simulations' accuracy we also performed the 
same measurements with some non-simulated faint guide stars. 

Figure |9l shows, at a 10% selection, the reference star flux 
requirements at two seeings. The requirements are similar at all 
percentage selections. 

At a guide star flux of 150 photons/frame (I — hi 6m) the 
resolution is reduced by 22% from that achieved with a very 
bright star. Because there appears to be a rapid falloff in reso- 
lution below that flux level we adopt 150 photons/frame as the 
minimum flux required for high resolution imaging. 

The more spread out PSFs in poorer seeing require more 
photons to give an acceptable SNR. The real (non flux rescaled) 
guide stars give matching results, although with higher noise 
induced by anisoplanatism. 

The ~150 photons/frame requirement can be understood in 
the following way. Even in poor seeing the best frames consist 
of a single speckle surrounded by a halo. If the speckle con- 
tains only a few tens of photons, Poisson shot-noise becomes a 
limiting factor for frame selection on the basis of Strehl ratios 
- thus degrading the output resolution. 

If a field contains several stars which are individually 
fainter than the guide star limit we have found it is possible 
to add the stars' PSFs together to provide a useful estimate of 
a bright guide star PSF, which can then be used in the standard 
way for Lucky Imaging. 

6.2. Sky Coverage 

The isoplanatic patch and reference star flux requirements di- 
rectly give us LuckyCam's achievable sky coverage. If the sys- 
tem is to achieve its most useful resolution gains we require an 
I ~ 16 reference star within 25" of our target. 

At a galactic latitude of 60 - 70° approximately 7% of the 
sky can be used with chance-placed guide stars, based on a ver- 
aged star counts from the USNO-Bl.O catalogue (Monet e t alJ 
2003). 30% of the sky is accessible at a latitude of 20 - 30° 
while at lower galactic latitudes sky coverage is virtually com- 
plete. 

6.3. Frame rates 

Freezing the effects of atmospheric turbulence requires expo- 
sures that are shorter than the atmospheric coherence time. One 
sufficiently fast (60Hz) trial dataset was taken 2003 and can be 
temporally rebinned to produce a range of effective effective 
frame rates (fig.llO>. Although this is a single observation the 
performance at 12FPS and 30FPS is consistent with that found 
in our earlier more numerous trials and is thus probably not 
atypical. 

Selection of frames (as opposed to simple shift-and-add) 
gives an improvement at all measured frame rates. At even 0.3 
FPS the frame selection gives gains in FWHM of over 1.6x, 
with an estimated limiting reference star magnitude of I=-i-19m 
in dark time. 

At ~30 FPS the atmospheric coherence time is well sam- 
pled in the best 10% of frames; there is no improvement with 
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Fig. 10. Imaging performance as a function of frame rate for in 
0.42" seeing. Shift-and-add (100% selection) FWHM perfor- 
mance is always a significant improvement over the autoguider 
performance because the autoguider uses the image centroid 
rather than the PSF peak as a position reference. 

faster frame rates. The frames with poorer quality do, however, 
benefit from still faster rates. The periods of poor seeing have 
a relatively smaller atmospheric coherence time because there 
are more isophase patches per unit area to be swept over the 
telescope aperture by the wind. The more inclusive frame se- 
lections include these periods and so require a faster frame rate 
than when selecting only the very best frames. 

Rimning at only 12FPS the current system cannot properly 
demonstrate the full potential of Lucky imaging; we have tested 
a new version of LuckyCam for science use in Summer 2005 
that provides the faster full-field frame rates required for full 
performance. 

7. Conclusions 

Lucky Imaging offers reliable and very valuable resolution en- 
hancement in all encountered conditions, even with our slow 
non-ideal camera. In poor seeing FWHMs are improved by up 
to a factor of four; in good seeing near diffraction limited im- 
ages are achieved. The angular radius over which the system 
can usefully operate in 1-band is as large or larger than our 
22"x20.5" field of view in all measured fields. The current sys- 
tem can work with guide stars as faint as I = 16 with full per- 
formance, giving near 100% sky coverage at galactic latitudes 
< 20° and >7% sky coverage at < 70°. Signal-to-noise ratio 
losses from frame selection are compensated by the greatly im- 
proved resolution. 

The system increases the range of conditions that are usable 
for particular observations - and can thus use otherwise lost 
telescope time, as well as giving resolutions in the visible far in 
excess of those of conventional CCD cameras on ground-based 
telescopes. 

Science programmes targetting both galactic (eg. binary 
candidates, brown dwarfs, globular cluster cores) and extra- 
galactic (eg. quasar host galaxies, damped Lyman-a absorbers) 



objects will benefit greatly from uncomplicated, reliable and 
cheap high resolution I-band imaging. 
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